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We have developed a new framework of Antisymmetrized Molecular Dynamics (AMD), to ade- 
quately treat the 7 = KN interaction, which is essential to study kaonic nuclei. The improved 
points are 1) pK~/nK° mixing and 2) total spin and isospin projections. These improvements enable 
us to investigate various kaonic nuclei (ppnK - , pppK - , pppnK - , 6 BeK - and 9 BK - ) systematically. 
We have found that they are deeply bound and extremely dense with a variety of shapes. 

PACS numbers: 



I. INTRODUCTION 

Recently, it was proposed that a K - meson can be 
deeply bound in light nuclei as a discrete state, such as 
3 He+K - , 4 He+K - and 8 Be+K - , where the K - meson 
makes the nucleus shrink drastically to form a dense state 
[J. Exotic proton-rich bound systems with K, such as 
ppK - , pppK - , pppnK - and 9 B+K - , are expected to be 
produced in (K - , 7r - ) reactions 0. In our previous pa- 
per, we investigated kaonic nuclei, which are denoted as 
K nuclei hereafter, 3 He+K - and 8 Be+K - , with a simple 
version of Antisymmetrized Molecular Dynamics (AMD) 
3]. Although our results are similar to those obtained 
in Ref. 0, a strange property appeared in 8 Be+K - . It 
is an isovector deformation, which means that the pro- 
ton distribution differs from the neutron one in spite of 
N = Z. Thus, K nuclei seem to provide interesting phe- 
nomena and stimulate further studies. 

Apparently, these interesting properties of K nuclei can 
be attributed to the KN interaction. Especially, the 1 = 
KN interaction plays an essential role. According to a 
precise experiment yj, the Is atomic state of a kaonic 
hydrogen (i.e. proton+K - ) is shown to be shifted up- 
ward. This upward shift, which is consistent with the 
low-energy scattering data of KN 5], suggests that the 
KN interaction is strongly attractive, so that the system 
of a K - and a proton has a nuclear bound state which 
corresponds to the 1 = A(1405) hyperon resonance ly- 
ing at 27 MeV below the K - p threshold. In a boson 
exchange potential model, the Julich group |(| showed 
that all of the cj, p and a mesons work coherently to 
give a strong attraction between a K and an N which ac- 
commodates a K - -p bound state, identified as A(1405). 
Studies based on chiral SU(3) Q also show that the 1 = 
KN interaction is attractive enough to form a A(1405). 
In this paper we employ a phenomenological KN interac- 
tion whose t-matrix in the KN channel is similar to 
that led by chiral SU(3) theory. 

Since the 1 = KN interaction is much more attrac- 
tive than the 1=1 one, a K - meson attracts protons 
rather than neutrons, causing an isovector deformation. 



Thus, the 1 = KN interaction is essential for studying 
K nuclei. 

In this paper we present systematic studies of K nu- 
clei with AMD. Since AMD treats a system in a fully 
microscopic way without any assumption concerning the 
structure of the system, it is suitable for studying K nu- 
clei, whose structures might be exotic. The simple ver- 
sion of AMD 3], however, has a technical problem in 
treating the 1 = KN interaction, which dominates K 
nuclear systems: it cannot include the degree of freedom 
of K°, and fails to describe the 1 = KN state. There- 
fore, in our previous paper we dealt with the 1 = KN 
interaction by incorporating its K°n part effectively into 
the K - p interaction. Of course, the 1 = KN interaction 
should be treated as precisely as possible, because it plays 
an essential role in K nuclei. In this paper, we improve 
the framework of AMD so that it can treat the 1 = 
KN interaction adequately. We introduce the degree of 
freedom of K° into the model space of AMD ("pK - /nK° 
mixing"). Since K nuclear states depend largely on their 
isospin due to the strong isospin-dependence of the KN 
interaction, we carry out the isospin projection as well as 
the angular momentum projection ("J & T projections") 
of the obtained intrinsic state. 

With the new version of AMD, we systematically in- 
vestigate a variety of K nuclei. We try to answer the 
following questions: i) What K nuclei are deeply bound 
with narrow widths? ii) Is there any strange structure 
peculiar to K nuclei? 

This paper is composed as follows: In Section |H] we 
present the improvements of AMD; pK - /nK° mixing, 
J & T projections, and other formalisms. In Section 
IIIII we demonstrate the capability of our new framework, 
and then apply it to various K nuclei (ppnK - , pppK - , 
pppnK - , 6 BeK - , and 9 BK - ). The results and a discus- 
sion are given in Section HVl We summarize our study in 
Section |V| 
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II. FORMALISM 

In the present study, we employ the AMD approach as 
a means of studying K nuclei. It has succeeded in study- 
ing the structures of light unstable nuclei |8| . In particu- 
lar, it is powerful for investigating the global properties of 
many light nuclei systematically. We know that various 
few-body methods, such as summarized in 9], can treat 
few-body K nuclei more accurately than AMD. Com- 
pared with these usual methods, AMD is more handy 
and applicable to more complex nuclei. Our aim is a 
systematic study of a variety of K nuclei. 

In the simple version of AMD employed in our previ- 
ous study, we restricted its model space to the proton, 
neutron and K - meson. Due to the lack of K° in the 
model space, we could not describe the 1 = KN state, 



\KN(I = 0)> 



1 

71 



(\pK-) + \nK }) . 



(1) 



In other words, we could not precisely treat the coupling 
of the pK~ pair with the nK° one through the 1 = 
KN interaction in the particle basis treatment of AMD. 
In the previous study, we incorporated all contributions 
from the 1 = KN interaction into an effective K _ p 
interaction, as follows: 



K-p 



a Vr 



KN(I=0) 



■pVj 



KN(I=1) •> 



(2) 



where a and /3 are some constants determined by count- 
ing the number of / = pairs and 1 = 1 ones in 
a given state of total isospin T. For example, we set 
(a,/?) = (§, j) in the case of ppnK - (T = 0). However, 
we have to check how this prescription is reliable for var- 
ious K nucleus cases. For this purpose, we introduce the 
degree of freedom of K° into the AMD framework to treat 
"pK _ /nK° mixing" directly. 



A. pK /nK° mixing 

First, we explain our idea for the simple case of a K 
nucleus ppnK - . In this K nucleus, the component of 
pnnK is mixed because a pair of pK - is replaced with 
that of nK° by the 1 = KN interaction. Hereafter, we 
express this state as |.H. An ordinary way to treat it is 
to perform a coupled channel calculation, preparing sev- 
eral Slater determinants for both channels of ppnK - and 
pnnK . In the present paper, we deal with such systems 
where several channels are coupled as follows: In stead 
of multi Slater determinants, we employ a single Slater 
determinant with charge-mixed single particle wave func- 
tions, i.e. 



\Ni) = Xi\p) + yi \n), 
\K) = z\K-) + w\K°), 



(3) 
(4) 



where \Ni) and \K) indicate a single nucleon wave func- 
tion and a K meson wave function, respectively. \Ni) 



can describe the state where a proton and a neutron are 
mixed, and also \K) can describe the state where K - and 
K° are mixed. With these wave functions we describe 
|.H as \det[N\N2N^]K). This state contains the com- 
ponent of \pnnK°) as well as that of \ppnK~). In this 
method, since each nucleon has a chance to be a proton 
or a neutron, the important channel is automatically cho- 
sen in the process of the energy variation. In addition, 
we expect that | det[NiN2Ns]K) can represent a state in 
which the contribution of several configurations is coher- 
ently additive: for example, if two configurations such 
as \ppnK~) and \pnnK°) work coherently, such state is 
represented as \p (p + n)n (K~ + K )). 

However, we remark one point: | det[NiN2Ns]K) is 
likely to have incorrect components, for example pppK - , 
ppnK and so on, which should not couple with |.H. 
To avoid the mixing of such incorrect components, we 
project it onto a state whose isospin-z component T z is 
equal to that of |.H. 

Now, we show the details of our wave function based 
on the concept of pK - /nK° mixing. Our nucleon wave 
function, \<pi), is represented by the superposition of sev- 
eral Gaussian wave packets 



l Dy 



fa) =2^ c '> x p 



WM). (5) 



Namely, the i-th nucleon is described by the superpo- 
sition of N n Gaussian wave packets whose centers are 
{Z % a }. \<Ji) means a spin wave function, and is | |) or 
||). \rl) means an isospin wave function, and has the 
following form: 



K) 



2 lc 



\P) 



(6) 



where 7^ is a variational parameter. In the usual AMD 
the isospin of each nucleon does not change, i.e. in the 
process of energy- variation the protons remain as protons 
and the neutrons as neutrons. However, in the present 
paper we make the isospins of all nucleons changeable so 
that we can treat pK~/nK° mixing. In the same way, a 
K meson wave function, \<Pk), nas the form 



N K 

\<Pk) =J2C« 



exp 



7K 



K). (7) 



Here, the isospin wave function of a K, \r^), as well as 
that of a nucleon, is changeable, 

|rf)=Q+7f)l^°)+Q-7f)l^-). (8) 

Because a nucleon is a fermion, we antisymmetrize 
the wave function of the nucleon's system, \$>n) — 
det[|(/^(j))] . Then, the K meson wave function is com- 
bined to it, |$) = I^at) \<Pk)- Moreover, we project the 
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total wave function onto the eigen-state of parity, 

I**) = -L[|$)± m ]. (9) 

If we perform an energy-variation with a trial wave 
function it is likely that the z-component of the 
isospin (T z ) of the total system is different from that of 
a K nucleus that we try to calculate originally. To avoid 
any mixing of components having an incorrect T z , we 
project the total system onto an eigen-state of T z before 
the energy-variation, 

\Pm^) = JdO exp[-iO(f z - M)]|$±). (10) 

Thus, we can obtain a wave function containing only the 
components of T z = M. We utilize |Pm < ^ )± ) as a trial 
wave function. 

Our wave function includes complex number parame- 
ters {X* a } = {(%, Zlri; C*, 7 *} and a real num- 
ber parameter (u). These are determined by the frictional 
cooling equation, as mentioned in S flll Dl 

In the present study, we use a common width param- 
eter (y) of a Gaussian wave packet for a nucleon and for 
a K meson so as to simplify our calculation. However, 
it seems natural that the spreading width of a nucleon 
is different from that of a K meson. We take this point 
into account by using different numbers of Gaussian wave 
packets for a nucleon and a K meson, i.e. N n in Eq. JSJ) 
is not equal to Nk in Eq. J7J). 

B. J &; T projections 

The angular momentum projection (J projection) is 
necessary to study K nuclei as well as that of usual nu- 
clei. In addition, the isospin projection (T projection) 
also seems to be important because the KN interaction 
has strong isospin-dependence. Therefore, we perform 
angular-momentum and isospin projections simultane- 
ously. J projection is done numerically by rotating the 
system in space, as has often been done. T projection 
is performed in quite the same way, but by rotating in 
isospin space. Our J & T projections are as follows: 

\^MK PtzT'^^ 1 ) = J dQang.DMK^ang.)RangX^ang.) 

x J dQ iso D^* T ,(n iso .)R iso X^iso.) 1$*), (11) 

where is the intrinsic wave function, which is al- 

ready determined by the frictional cooling equation, as 
shown in S flll Dl We calculate various expectation values 
with \P J MK P£ T , 

C. Hamiltonian 

Our Hamiltonian in AMD calculations, 

H = T + V NN +V c + V K N-f G , (12) 



is composed of the kinetic energy T, the effective NN 
potential Vat at, the Coulomb force Vc, and the effective 
KN potential Vkn- The center-of-mass motion energy, 
Tg, is subtracted. In the kinetic energy and the center-of- 
mass motion energy, we treat the mass difference between 
a nucleon and a K meson correctly. For example, the 
kinetic energy is 

f=y^ + ^, (is) 

1=1 

where tun and tuk indicate the mass of a nucleon and 
that of K~, respectively. The Coulomb force is rep- 
resented by the superposition of seven-range Gaussians 

In the study of K nuclei, we do not use existing ef- 
fective interactions which may be justified for studying 
phenomena around the normal density. Since the system 
is likely to become extremely dense due to the strong 
K attraction, we employ the ^-matrix method 0. We 
adopt the Tamagaki potential (OPEG) 12] as a bare NN 
interaction, and the Akaishi-Yamazaki KN potential Q 
as a bare KN interaction. Because the Tamagaki poten- 
tial can reproduce NN phase shifts up to 660 MeV in the 
lab. 12], we expect that it can be applied to such ex- 
tremely dense states. The effective NN/KN interactions 
constructed from the bare ones are represented by the 
following ten-range Gaussians: 

10 

V$ N {r) = Y, V NN, a eM-(r/r a ) 2 } (14) 

a=l 
10 

VR N (r) = Y, V RN, a eM-(r/r a ) 2 } (15) 

a=l 

"X" in Eq. jnj) is 3 E, 1 or 3 0, and "F in Eq. (fTKl) 

is or 1. We use these V^ N (r) and Vj^ N (r) as effective 
NN/KN interactions in our AMD calculation. 

Our procedure is as follows: 1) For a given density 
and starting energy of K~, we construct a ^-matrix. 2) 
Using the ^-matrix we carry out the AMD calculation. 3) 
After the AMD calculation, we check whether or not the 
obtained density and binding energy of K~ are consistent 
with those of the ^-matrix used in the calculation. 4) If 
no consistency is accomplished, we guess and impose a 
new density and a new starting energy of K~ for the g- 
matrix calculation by referring to the results obtained so 
far and return to 1). We repeat this cycle until obtaining 
a consistent result. 



D. Frictional cooling equation with constraint 

Our wave function contains complex variational pa- 

rametera, {X*} = {C£,Zj,, 7 i ; Cf,z£ 7 f}. They 
are determined by the energy-variation. In our study, 
we employ the frictional cooling method as a means of 
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TABLE I: Results of ppnK - in various N n and N K . B.E.: 
total binding energy. p(0): central density. Rrms, Rrms, 
Rrms'- root-mean-square radii of matter, proton and neutron, 
respectively. (/3, 7): deformation parameters. 



N n 


N K 


B.E. 


P(0) 


Rrms 


TDP 

J L rms 


-Ll rmS 


P 
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[MeV] 


[fm- 3 ] 


[fm] 


[fm] 


[fm] 




[deg.] 


2 


5 


112.7 


1.39 


0.72 


0.70 


0.75 


0.19 


0.0 


3 


5 


113.8 


1.40 


0.73 


0.70 


0.77 


0.14 


49.0 


4 


5 


117.0 


1.41 


0.72 


0.69 


0.78 


0.02 


49.4 


2 


8 


113.5 


1.37 


0.72 


0.71 


0.75 


0.19 


0.0 


2 


10 


113.7 


1.37 


0.72 


0.71 


0.75 


0.18 


0.0 


3 


10 


115.1 


1.41 


0.72 


0.70 


0.77 


0.14 


51.7 


4 


10 


116.9 


1.49 


0.71 


0.68 


0.77 


0.00 


54.4 



energy- variation, 



dH 

dxj: 



rj 



and C.C.(16) 



Here, TL is the expectation value of the Hamiltonian and 
W is a constraint condition. 77 is a Lagrange multiplier, 
which is determined by dW/dt = 0. It is easily proved 
that, if we assume /i < in Eq. (fTrlj) and all of the param- 
eters are developed with time according to Eq. (p"6|) . the 
energy of the system decreases while satisfying the con- 
straint condition W = 0. If we use the superposition of 
several Gaussian wave packets to represent a nucleon and 
a K meson wave function, we need a constraint condition 
in order to fix the center of mass of the total system to 
the origin, and then W is expressed as follows: 



(R G f + {P G f 



Rg 



Yh=i m NTi + m K r K 



A m at + rriK 



G 



Pk- 



(17) 
(18) 

(19) 



III. TESTS OF OUR METHOD 

Before applying our method to studies of various K 
nuclei, we investigate the basic properties of our method. 



we find that the total binding energy and the central den- 
sity are almost converged up to N n — 4 and Nk = 10. 
However, we notice that the shape of the system, repre- 
sented by the deformation parameters (/?, 7), is strongly 
dependent of N n . This phenomenon can be understood 
as follows. As mentioned in our previous study 3], the 
protons distribute compactly near a K - so as to decrease 
their total energy by the strongly attractive K - p inter- 
action. On the other hand, the neutron is widely spread 
and its total energy decreases by reducing its kinetic en- 
ergy. Therefore, the protons stay compactly inside the 
system, while the neutron remains widely outside of the 
system. Thus, the neutron contributes to the shape of the 
total system. In the case of N n = 2, since the neutron is 
represented by two Gaussian wave packets, it can spread 
only linearly. Thus, the total system deforms prolately. 
In the case of N n = 3, it can spread with a triangular 
shape. Thus, the total system deforms oblately. In the 
case of N n = 4, it can spread with a tetrahedron shape. 
The total system is therefore spherical. Thus, the shape 
of the total system changes as N n is varied. 

Such a dependence of the shape on 7V n seems to be 
peculiar to ppnK - where the neutron number is equal 
to 1 and the proton number is 2. In addition, the to- 
tal binding energy and the central density do not so 
strongly depend on N n and Nk- Therefore, taking the 
cost-performance of calculations into account, we adopt 
the model space of N n = 2 and Nk = 5 in our calcula- 
tions. 



B. Solution of ppnK 

We now check whether our new framework, pK - /nK° 
mixing and J & T projections, works correctly or not. 
We perform a test on a system of ppnK - . 

First, we investigate the property of J & T projections. 
Although only the J projection has often been carried out 
in the study of light unstable nuclei 0, , the present 
study for the first time makes the T projection. In Ta- 
ble II, we show various quantum numbers of the wave 
function before projection (|Pm < ^ )± )) and that after pro- 
jection (\Pm K P^ T /^ >± )). Apparently, the ground state 

of ppnK - seems to have quantum numbers of J n = ^ + 
and T = 0. We performed J & T projections so that 
the total system had such quantum numbers. Table II 



A. Dependence on the number of wave packets 

As shown in Eq. JSJ) and Eq. J7J) , we represent a single 
nucleon wave function and a K meson wave function with 
N n and Nk Gaussian wave packets, respectively. We in- 
vestigate how much the solution depends on N n and Nk- 
First, we perform a test in the case of ppnK - without 
pK - /nK° mixing for simplicity. Table I shows the re- 
sults of ppnK - for various Af n and Nk- From this table, 



TABLE II: Quantum numbers before and after projection, 
jj: (J ) of total system. J^, L 2 N and S%: (J ), (L ), and 
(S ) of nucleon system. L 2 K \ {L ) of a K meson. T 2 and T z : 
(f 2 ) and if 2 ) of total system. 







T 2 


T 2 




T 2 


rj~\2 


T z 


After 


0.75 


0.78 


0.03 


0.75 


0.03 


0.00 


0.00 


Before 


1.36 


1.22 


0.44 


0.78 


0.14 


0.02 


0.00 
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TABLE III: Results of ppnK - . B.E.: binding energy mea- 
sured from the threshold of ppn+K - . F: width decaying to 
Aty. p(0): central density. R r ms' root-mean-square radius of 
nucleon system. 





B.E. 


r 


P(0) 


Rrms 




[MeV] 


[MeV] 


[fm -3 ] 


[fm] 


Present 


110.3 


21.2 


1.50 


0.72 


simple AMD 


105.2 


23.7 


1.39 


0.72 


BHF 


108 


20 




0.97 



cle base: i) V nK - and V^o, ii) V pK - and F n ^o, and iii) 
VpK- ,nK°- The interactions i) and ii) are working in each 
channel of ppnK - and pnnK , and their expectation val- 
ues are equal to —45 MeV and —255 MeV, respectively. 
Interaction iii) is related to pK - /nK° mixing through 
the 1 = KN interaction, and its expectation value is 
equal to —88 MeV. Thus, we find that the binding of 
this system is mainly due to the type ii) interaction and 
is further supported by the type iii) interaction, which 
causes coupling between the two channels. 



shows that (J ) = 0.75 and (f ) = 0.00, which agree 
with J( J + 1) = \ • § and T(T + 1) = • 1, respectively. 
Therefore, it is found that our J & T projections work 
well. 

Next, in Table III various quantities obtained in the 
present calculation (Present) are compared with our pre- 
vious result of a simple version of AMD (simple AMD) 
and the result of a BHF calculation (BHF) 1]. This 
table shows that the present result is almost identical to 
others. Since the isospin-z component of each particle 
is changeable in the present framework, we investigated 
each particle-number. Although we calculated ppnK - , 
the numbers of protons and neutrons after the calcula- 
tion are both equal to 1.5, while those of K - and K° are 
both 0.5. This means that ppnK - and pnnK are mixed 
with a ratio of 1:1 as the result of pK - /nK° coupling 
through the 1 = KN interaction. 

Here, we remark on the components of the KN inter- 
action. We can separate it into three parts in the parti- 



C. Interpretation of the density distribution 

The influence of pK - /nK° mixing can clearly be seen 
in the density distribution. Fig. [I] displays the density 
distribution of protons and neutrons in the ppnK - sys- 
tem calculated by the new framework. We can see that 
the proton distribution is almost the same as the neutron 
one, contrary to our previous result Q, where protons 
distribute more compactly than neutrons because of the 
strong attraction between the K - and the proton. We 
can solve this contradiction by introducing the concept 
of an intrinsic state in isospin space, as follows: the cal- 
culated expectation value of T 2 with the state \Pt z =o 
in the case of ppnK - , is nearly equal to zero. Therefore, 
this state is the eigen-state of isospin, i.e. T = 0, and we 
express it as |^H (T = 0)) hereafter. It is easily found 
that ||-H (T = 0)) is composed of two configurations 
concerning the z-component of isospin: 



||H(T = 0)> 



Pt z =o [\&n) <8> \vk}} 

%= 



Pt z =o 

P t n = i $AT 



m= — oo 



T K._ 



$n) ® ( Y ^ T Z K -- 

m=±l/2 



Vk) 



(20) 



where P t n and P t k are TVprojection operators for the 
nucleon system and the K meson, respectively. Accord- 
ing to the values of and , the first term indi- 
cates ppnK - , while the second term indicates pnnK . 
Hereafter, we express them as \ppnK~) and \pnnK°), 
respectively. In addition, the overlap between \pnnK°) 
and e l7vTy \ppnK~), calculated numerically, is found to 
be almost one. Since this fact indicates \pnnK°) c± 
e tnTy \ppnK~), we can say that our wave function sat- 
isfies 

i|h(t = o)>= J2 etefy \pp nK ~) 

6> = 0,7T 



= e%dfv \vnnK*). (21) 

(9 = 0, 7T 

Namely, the \ppnK~) state rotates in isospin space so 
that ||-H(T = 0)), which has the good quantum number, 
T = 0, is formed. Based on the analogy of an "intrinsic" 
state in a deformed nucleus, rotating in the space to form 
the eigen-state of the angular momentum, we can regard 
the \ppnK~) or equivalent \pnnK°) as an intrinsic state 
in the isospin space of ||-H(T = 0)). 

Fig. [21 displays the proton and neutron distributions 
of the intrinsic state \ppnK~). Clearly, the proton distri- 
bution is more compact than the neutron one; this fact 
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FIG. 1: Proton (left) and neutron (right) distributions of 
ppnK - obtained in the present calculation. The size of each 
frame is 3fm x 3fm. 
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FIG. 2: Proton (left) and neutron (right) distributions of 
ppnK - . See the text for details. 



is consistent with our previous study. Here, we note one 
point. Investigating the density distribution of \pnnK°), 
we find that its proton and neutron distributions are com- 
pletely identical with the neutron and proton ones in the 
\ppnK~), respectively. In other words, neutron distribu- 
tion is more compact than proton one in the \pnnK°). 
After all, since Fig. [I] is drawn with the ||-H(T = 0)), 
which includes the \ppnK~) and the \pnnK°) to form the 
eigen-state of isospin, the proton distribution is quite the 
same as the neutron one. 



IV. RESULTS 

Now that we have confirmed that the new framework 
works well, we proceed to an investigation of pppK - , 
pppnK - , 6 BeK - and 9 BK - with this framework. Since 
the intrinsic state of isospin space is found to be mean- 
ingful, as mentioned in the previous section, we generally 
denote K nuclei by their intrinsic states without any loss 
of validity; for example, ~.H is represented by ppnK - . 
In all calculations, a single nucleon and K meson are 
described with two Gaussian wave packets and five, re- 
spectively, namely N n = 2 in Eq. © and Nk = 5 in Eq. 

Q 



TABLE IV: Summary of present calculations. J 7 * and T: 
spin-parity and isospin of total system. E(K): binding en- 
ergy of K - meson. T: width decaying to £7r and Aty chan- 
nels. p(0): central density. R r ms- root- mean- square radius of 
nucleon system. (/3, 7): deformation parameters. 





r 


T 


E(K) 


r 


P(0) 


Rrms 


P 
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[MeV] 


[MeV] 


[fm -3 ] 


[fm] 




[deg.] 


ppnK 


1 + 
2 





110.3 


21.2 


1.50 


0.72 


0.22 


9.2 


PPPK" 


3 - 
2 


1 


96.7 


12.5 


1.56 


0.81 


0.70 


11.8 


pppnK - 


1~ 


1 
2 


105.0 


25.9 


1.29 


0.97 


0.54 


3.8 


6 BeK - 


0+ 


1 
2 


104.2 


33.3 


0.91 


1.17 


0.44 


0.3 


9 BK- 


3 - 
2 





118.5 


33.0 


0.71 


1.45 


0.46 


20.8 



we confirm after J & T projections that the expectation 
values of J and T are equal to those we set beforehand. 
Fig. [31 displays the behavior of the total binding energy 
and the decay width, from ppnK - to 9 BK - . The previ- 
ous result 0| of 8 BeK - is also shown there. According to 
the Table IV, all K nuclei are bound by about 100 MeV. 
Fig. 01 shows us that the K nuclei except for pppK - are 
bound below the £7r threshold which is the main decay 
channel. 

In Fig. the thresholds of various decay modes are 
also shown, from which we can estimate the stability 
of pppnK - to 9 BK - for the strong interaction. First, 
pppnK - is found to be unstable for the strong interac- 
tion, since its total binding energy (112.5 MeV) is smaller 
than that of ppnK - + p (117.8 MeV). On the other hand, 
the deepest thresholds for 6 BeK - , 8 BeK - and 9 BK - are 
ppnK-+ 3 He (125.3 MeV), ppnK - +a+n (143.5 MeV) 
and 8 BeK - +p (159.0 MeV), while their total binding en- 
ergies are 131.7 MeV, 159.0 MeV and 173.5 MeV, respec- 
tively. Therefore, 6 BeK - , 8 BeK - and 9 BK - are stable 
for the strong interaction. 

As mentioned above, the K - meson is bound by about 
100 MeV in all of the K nuclei that we calculated. pK - 
is bound by 27 MeV based on the initial assertion that it 
forms A(1405), and ppK - is bound by 48 MeV according 
to |2|. Therefore, the binding energy of the K - meson 
seems to be saturated in K nuclei heavier than ppnK - . 
We think that this saturation of the binding energy is 
related to the range of the KN interaction. Since it is 
very short, the number of nucleons which a single K - 
meson can interact with is limited. We count up the nu- 
cleons staying in the region where the K meson's density 
falls down from the maximum value Pmax t° \ Pmax- 



A. Binding energies 

The results are summarized in Table IV. We determine 
J n and T by assuming that nucleons occupy the configu- 
ration of the normal ground state and a K meson occupies 
the 0s state. Of course, as done in the former section, 



TABLE V: Binding energy and number of strongly interacting 
nucleons near the K meson. 





ppnK 


pppK" 


pppnK 


6 BeK~ 


9 BK- 


B.E. 


110.3 


96.7 


105.0 


104.2 


118.5 


Nucleon 


1.67 


1.14 


1.78 


2.55 


2.53 
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FIG. 3: Behavior of the total binding energy and the decay width, from ppnK - to 9 BK - . The decay width is to Att and E7r 
channels. The red point and line indicates the binding energy and width of the K nucleus. The blue-dashed line corresponds 
to the S7T threshold. The thresholds for other decay modes are expressed by the green-dashed lines. 
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FIG. 5: Proton (left) and neutron (right) distributions of 
6 BeK~. 

In Table V, we show the number of strongly interacting 
nucleons ( "Nucleon" ) staying around the K in various K 
nuclei. Except for pppK - , which has a peculiar struc- 
ture, as mentioned in the later section, about 1.7 to 2.5 
nucleons are found to stay near the K meson. 

B. Density distribution 

Fig. 0] displays the nucleon density distributions of 
ppnK - , pppK - , pppnK - , 6 BeK - and 9 BK - . It is found 
that K nuclei have extremely dense and peculiar nucleon 
distributions. 

6 BeK - has a two-center-like structure similar to 



8 BeK - Q. Fig. [5] shows proton and neutron distribu- 
tions separately. We can find that protons have a two- 
center-like structure and that neutrons stay between two 
pairs of protons, against our expectation that 6 Be should 
have such a structure as a+2p. We note that the K me- 
son's density distribution, which is not displayed here, is 
very similar to the neutrons' one. The structure of 9 BK - 
is quite similar to that of 8 BeK - |3|. 

The most exotic system of pppK - shows a very pecu- 
liar density distribution. Strictly speaking, this system 
has not only the component of pppK - , but also that of 
ppnK due to the 1 = KN interaction. We express this 
system as |-He. Fig. shows only proton density distri- 
bution extracted. We can see a "satellite" in this figure. 
Summing up the density of proton in the region of this 
satellite, we find that the proton number is nearly equal 
to one. Thus, this satellite is a single proton. In addi- 
tion, the particle numbers of proton, neutron, K - and K° 
are 2.67, 0.33, 0.67 and 0.33, respectively. We can under- 
stand these particle numbers and the density distribution 
consistently as follows. We regard this system as a single 
proton combining |.H, following its density distribution 
(Fig. |SJ). In isospin space, |.H means the \T = ~, T z = |) 
state, which is composed of two nucleons and a K meson. 
The weight of each component included in it is deter- 
mined by Clebsch-Gordan coefficients as shown below: 



||He) = |p)®&H> 
= \p) <g> 

= \p) ® I \l i \pp ® K 



T—-T - - 

" 2' 2 2 



T N = 1 T N = 1T K = i T K 



\pn K°) 



T N = l T N _ q jiK _ 1 T K = 1 

' z ' 2 z 2 



(22) 



The particle numbers counted according to Eq. ((221) are 
quite the same as the above values. After all, we come 
to an idea that in |-He one proton keeps its identity and 
that the residual part is composed of ppK~ and pnK°, 
which are mixed due to the 1 = KN interaction. Note, 
however, that the single proton is strongly bound to |-H. 
In addition, we studied the dependence of this system on 
the number of wave packets. Even if one nucleon is repre- 
sented by four Gaussian wave packets (i.e. 7V n = 4 in Eq. 
(J5J)), a proton is still pushed out, but a little less clearly. 
We consider that the third proton is pushed up to the Op- 
shell due to Pauli Blocking, so that it forms a satellite- 
like structure. pppnK - has an extra neutron compared 
to pppK - . When a neutron is added to pppK - , such a 
satellite-like structure disappears and a different struc- 
ture is formed for pppnK - . 



V. SUMMARY 

We improved the antisymmetrized molecular dynam- 
ics (AMD) regarding two points and applied it to sys- 
tematic studies of K nuclei. One of our improvements is 
"pK _ /nK° mixing", which enables us to treat directly 
the coupling of pK _ and nK° through the 1 = KN 
interaction. The other one is "J & T projections" with 
which the strong isospin dependence of the KN interac- 
tion is expected to be correctly treated. Thus, AMD is 
capable of treating the 1 = KN interaction adequately, 
which plays an essential role in K nuclei. 

We have investigated the properties of our new frame- 
work on ppnK - . After J & T projections, the total sys- 
tem is found to possess the quantum numbers J and T, 
which turn out to be equal to those which we had as- 
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FIG. 6: Proton distribution in pppK . 

signed beforehand. Namely, we have confirmed that the 
J & T projections work correctly. The new result and the 
previous one are very similar to each other, but we can 
see the influence of pK~/nK° mixing in the density dis- 
tributions of the protons and neutrons. Owing to the in- 
troduction of K°, the present wave function can form the 
eigen-state of isospin, i.e. ||-H(T = 0)). When we draw 
the proton and neutron distributions of ||-H(T = 0)), 
they are clearly different from those shown in our previ- 
ous study. We confirm that the ||-H(T = 0)) is formed 
by the rotation of |ppnK~) in isospin space. Namely, 
|ppnK~) is considered to be an "intrinsic state in isospin 
space" for the ||-H(T = 0)). This intrinsic state |ppnK~) 
is found to have quite the same proton and neutron distri- 
butions as those in our previous study. The ||-H(T = 0)) 

contains both |pnnK°) and |ppnK~) components with 
the same ratio. This is consistent with the fact that the 
calculated proton and neutron numbers are both equal 
to 1.5. The coupling of these two components due to 
the 1 = KN interaction helps the binding of the total 
system ||H(T = 0)). 

We have studied ppnK - , pppK - , pppnK - , 6 BeK - 
and 9 BK - with our new framework. All K nuclei that 
we investigated are bound by about 100 MeV below the 
threshold of each nucleus+K - . Except for pppK - , they 
are bound below the E-zr threshold, which is the main 
decay channel. Since their decay width is 20 to 40 MeV 



and small compared to their binding energy, they ap- 
pear to be discrete states. For the strong interaction, 
pppnK - is found to be unstable, while 6 BeK - , 8 BeK - 
and 9 BK - are stable. We found that they have very 
different structures. Especially, pppK - shows an inter- 
esting satellite-like which is composed of a single proton. 
This proton keeps its identity and is strongly bound by 
the main body. 

According to our present study, we predict various 
deeply bound K nuclei. They have very peculiar struc- 
tures with extremely high densities, which we have never 
seen. In the future, such K nuclei as those investigated in 
the present paper may be explored experimentally. For 
instance, ppnK - can be formed from a 4 He(stopped K - , 
n) experiment, which is now under way at KEK [JEH. 
The use of in-flight (K - ,Af) reactions is also proposed 
. Proton-rich exotic K nuclei are expected to be pro- 
duced by (K - , 7T~) reactions via A* doorway states [3- 

Finally, we would like to mention possibility of the new 
framework of AMD. Since it can be extended straight- 
forwardly to the case of multi K's, we can investigate 
multi-K nuclei, which are closely related to kaon conden- 
sation and strange quark matter. The success of the new 
framework of AMD means that we can deal with even 
more fields of physics, because coupled-channel-like cal- 
culation can be carried out in the new version of AMD. 
For example, we can study A-E mixing in hyper nuclei 
[l5| with this new framework. 
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